Revisiting the valence-band and core-level photoemission spectra of NiO 
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We have re-examined the valence-band (VB) and core-level electronic structure of NiO by means 
of hard and soft x-ray photoemission spectroscopies. The spectral weight of the lowest energy state 
was found to be enhanced in the bulk sensitive Ni 2p core-level spectrum. A configuration-interaction 
model including a bound state screening has shown agreements with the core-level spectrum, and off- 
and on-resonance VB spectra. These results identify the lowest energy states in the core-level and 
VB spectra as the Zhang-Rice (ZR) doublet bound states, consistent with the spin-fermion model 
and recent ab initio calculations within dynamical mean- field theory (LDA+DMFT). The results 
indicate that the ZR character first ionization (the lowest hole-addition) states are responsible for 
transport properties in NiO and doped NiO. 

PACS numbers: 71.10.-w, 71.28.+d, 79.60.-i 



NiO has played a very important role in clarifying the 
electronic structure of transition metal (TM) oxides. It 
has set a framework for understanding the rich physi- 
cal properties of TM compounds. More recently, the 
discoveries of (i) bistable resistance switching [1] (ii) a 
giant low frequency dielectric constant [2] and (iii) the 
direct correlation between exchange bias and interfacial 
spins have led to a major resurgence of interest in NiO, 
as well as doped and deficient NiO. In order to eluci- 
date the transport properties of NiO and doped or de- 
ficient NiO, one must understand the character of holes 
that enter the acceptor levels or the top of the valence 
band (VB). For a long time, NiO was thought to be a 
prototype Mott insulator in which the insulating gap is 
caused by on-site Coulomb energy U{= Udd)- Zaanen, 
Sawatzky and Allen redefined NiO as a charge- transfer 
(CT) insulator [4]. New theoretical concepts coupled with 
photoemission spectroscopy (PES)[5, 0, S, played a 
pivotal role in setting this stage. Motivated by resonant 
VB PESQ, Fujimori and Minami originally applied a 
configuration-interaction (CI) model to NiOfi]. It was 
shown that the first ionization state is of mainly 3(i^L 
character and Sd^-like states occur at a higher binding 
energy [6]. This picture indicated that the O 2p band ap- 
pears between the lower and upper Hubbard bands, and 
the insulating gap is formed between the O 2p and the 
upper Hubbard d bands. This approach is well estab- 
lished for explaining core-level PES of TM compounds in 
general. However, for NiO itself, the double peak struc- 
ture of the lowest energy states in the Ni 2p core-level 
PES cannot be explained by a single-site NiOe CI model. 
This shortcoming was overcome by using a multisite clus- 
ter (MSC) model[10], which includes nonlocal screening 
effects from neighboring Ni sites. It was shown that an 
extra peak appears at the high binding energy side of 



the main line due to the nonlocal screening. Because 
the structures seen in core-level spectra arise from ener- 
getically different screening channels in final states, this 
extra intersite screening channel should also be impor- 
tant for other spectroscopies, and particularly so for the 
VB PES. In fact, the similarity in the spectral shapes of 
the VB and the core-level is well-known [9]. Nonethe- 
less, the importance of the nonlocal screening in VB 
PES of NiO based on a MSC model has not been ex- 
plored to date. Theoretical studies of the VB electronic 
structure have shown that the two lowest energy elec- 
tron removal states are Zhang- Rice (ZR) doublet bound 



states pjj, ll2|, ll3|, [IJ, ll5|. They arise from a competi- 



tion between O 2j9-Ni 3(i hybridization and the Ni on-site 
Coulomb interaction. These results are consistent with 
angle resolved PES of NiO, which has shown nondisper- 
sive features for two lowest energy states [8]. While, for 
high Tc cuprates, the relation between the VB and core- 
level PES has been studied extensively in the context of 
the ZR singlet bound state, it has received less attention 
in NiO. This is partly because VB PES directly probes 
the electron removal spectral function {dP'^dP''^) while 
the core-level PES indirectly probes many body effects 
in the presence of a core-hole (2p^3(i^^2p^3(i^). This 
makes calculations much more difficult for NiO (d^^dJ), 
in the presence of multiplets and many body effects, as 
compared to the cuprates {d^^d^). Therefore, an ex- 
planation of the spectral shapes and its relation to the 
ZR doublet bound states in both the core-level and the 
VB, using a unified model, remains an open question. 
The main purpose of the present letter is to propose 
an explanation for both core-level and VB PES based 
on the CI model, in which the charge transfer from the 
bound state is included within a single-impurity frame- 
work. The model retains the simplicity and generality of 
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FIG. 1: (Color online) (a) Comparison between experimental 
Ni 2p core- level HX and SX PES. (b) A CI calculation for Ni 
2p spectrum. Final- state components are also displayed. The 
core-hole lifetime broadening is taken as 0.2 eV. The Gaussian 
broadening is 0.2 eV (half width at half maximum). 



the single cluster approach, while taking into account the 
screening effects due to the bound states in an approxi- 
mate way. 

Hard x-ray (HX) PES is a bulk sensitive probe of the 
electronic structure due to its ability to overcome high- 
surface sensitivity of conventional PES. Unlike soft x- 
ray (SX) PES, TM 2p core-level HX PES have shown 
additional well-screened features with significant inten- 
sity at the low binding energy side of the main peak[16]. 
It should be noted that these spectral changes are seen 
only in the metallic phase, and the features are weak or 
missing in the insulating phase. These features were ex- 
plained well by the CI model including a screening chan- 
nel derived from coherent states at Fermi energy. Recent 
theoretical studies of the MSG model 17, offer an al- 
ternative interpretation on the basis of a nonlocal screen- 
ing scenario. The question whether such features refiect 
the screening from the coherent states, or the nonlocal 
screening of the MSG model, and whether these models 
can be effectively mapped into each other or not, remains 
to be answered in general. However, at least in high Tc 
cuprates, the MSG model has confirmed that the screen- 
ing from the ZR singlet bound state at the top of the VB 
is the origin of the lower binding energy feature in the Gu 
2p core- level HX PES [17J . Because NiO is a good insula- 
tor, a big change in HX PES is not expected. However, 
the nonlocal screening effect was originally proposed for 
Ni 2p core-level PES of NiO. We, therefore, felt it worth- 
while to study it by HX PES. This is another motivation 
of the present letter. 

HX PES measurements were performed in a vacuum of 



1 X 10-^^ Torr at undulator beam line BL29XU, SPring- 
8 using a Scienta R4000- 1 Ok V electron analyzer. The 
instrumentation details are described in Ref. [3]. The 
energy width of incident x-ray was 70 meV, and the to- 
tal energy resolution, AE was set to ~ 0.2 eV. SX-PES 
and X-ray absorption spectroscopy (XAS) (total electron 
yield mode) were performed at BL17SU, with AE ~ 0.2 
eV. All PES measurements used normal emission geom- 
etry to maximize depth sensitivity. Sample temperature 
was controlled to ±2K during measurements. A single 
crystal of NiO was fractured in-situ. HX (SX) PES were 
measured at 200 K (room temperature) to avoid charg- 
ing. The Fermi level of gold was measured to calibrate 
the energy scale. The O Is spectra were clean single 
peaks, as is known for cleaved NiO. Gore- level spectra 
were corrected by subtracting a Shirley-type background. 

First, we present experimental Ni 2p core- level HX 
PES and SX-PES of NiO obtained with photon energies 
of hv=im keV and 1.25 keV in Fig. 1(a). The kinetic 
energy of the Ni 2p core-level corresponds to a probing 
depth of ~60 A for HX PES as determined by the in- 
elastic mean free path, while it corresponds to ~15 A for 
SX-PES [20]. The core- level spectra were normalized at 
the feature G, since it is well known that the feature G 
is insensitive to the surface l2l|, |22|. In comparison with 



SX-PES, the intensity of the feature A for bulk sensi- 
tive HX PES was found to be enhanced by a factor of 
1.35 with respect to the SX-PES (see inset). In contrast, 
negligible change was seen in the features B and G. In 
the MSG model, the most important result is that the 
feature B originates from the nonlocal screening. It was 
also predicted that nonlocal screening in the peak B is 
more important for the surface than for the bulk. How- 
ever, our experiments in Fig. 1(a) clearly show a different 
behavior: only the feature A was appreciably changed be- 
tween the HX and SX PES. It should be noted that the 
present SX-PES result shows the same qualitative be- 
havior as previous SX-PES results^, [i^], in which the 
intensity of feature A is significantly reduced at nearly 
grazing take-off angle i.e. in surface sensitive geometry. 
This clear enhancement of the feature A is therefore a 
direct consequence of the bulk sensitivity of HX PES. 

In order to address the origin of the observed behavior, 
we used a GI model with full multiplets in the Oh local 
symmetry and include a bound state to describe the 
spectra. The basic framework of our model is as follows: 
(i) we have used a single-site NiOe cluster with the usual 
screening channel for the charge transfer between the Ni 
site and the O sites, (ii) based on results of the spin- 
fermion model [lH and the LDA+DMFT calculation 
we have introduced charge transfer between the Ni site 
and the bound state within the GI model. The ground 
state is described by a linear combination of 3(i^, 3(i^I/, 
W^^I?, 3d^Z, Sd^^LZ, and Sd^^Z^, where Z and L 
represent a hole in the bound state and a hole in O 2p 
ligand states, respectively. The Hamiltonian is given 
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FIG. 2: (Golor online) (a) Experimental and calculated VB 
spectra with photon energies /izy=1.25 keV and 846 eV. (b) 
Final state components for the VB spectrum, (c) Experimen- 
tal and calculated resonant PES with (thick solid line) and 
without (dashed line) the bound state screening. Excitation 
energies are labeled on the L3 edge XAS (inset). 



by H=Hj-\-Hjj. The Hj represents the standard GI 

model, Hi=^^^esd{^)d'l^dra + E 

Era ^p(r)4a^ra+Era ^(r)(4a^ra+4a^ra)+i^muit 

- plna'P^(^')' Here SsdO^), S2p and ep{T) represent 
the energies of Ni Sd, Ni 2p and O 2p states, respec- 
tively, with the irreducible representation {T=eg and 
t2g) of the Oh symmetry. The indices m and a are 
the orbital and the spin states. V{T)^ [/, and —Udc 
are the hybridization between Ni 3d and O 2p states, 
the on-site repulsive Goulomb interaction between Ni 
3d states and the attractive 2p core-hole potential, 
respectively. The Hamiltonian i^muit describes the 
intra-atomic multiplet coupling between Ni 3d states 
and that between Ni 3d and Ni 2p states. The spin-orbit 
interactions for Ni 2p and 3d states are also included 
in i^muit- In addition to the usual GI model {Hj 
term), we have introduced the bound state labeled 'Z' 
as a new screening channel described by Hu term, 

Hn = Er.ec(r)4^cr. + Er.^*(r)(4.cra + 4,dr.). 
An effective coupling parameter V*(r) for describing 
the interaction strength between Ni 3d orbitals and the 
bound state is introduced, analogous to the hybridization 
y(r). The GT energy from the bound state to the Ni 
3d orbital is A*, whereas the usual GT energy A (from 
the O 2p states to the Ni 3d states) is defined as the 



energy difference of the configuration-averaged energies 
E{3d^L)—E{3d^). The parameter values (in eV) are 
summarized in Fig. 1. The Slater integrals and the spin 
orbit coupling constants are calculated by Hartree-Fock 
code [24] and the Slater integrals are scaled down to 85%. 

The calculated spectrum is shown in Fig. 1(b). The 
experimental features are reproduced well by the calcu- 
lation over the whole energy range. To clarify the peak 
assignment, we show in the lower panel of Fig. 1 the vari- 
ous final state components of the core-level spectra. It is 
clear that the dominant spectral weight in the features A 
arises from the 2p^3d^Z_ state, while the contribution of 
the 2p^3d^L is dominant for the feature B. The satellite 
G is mainly due to the 2p^3d^ state. The present as- 
signment is different from that by the MSG model [13] in 
which the peak A arises from the 2p^3d^L and the peak 
B is due to the nonlocal screening. 

In order to obtain a better understanding of the bound 
state screening, we applied the present model to VB 
PES. Figure 2(a) shows the calculated VB spectrum com- 
pared to the measured spectra with hv=1.2b keV and 
846 eV. Our 1.25 keV spectrum is consistent with previ- 
ous reports [1, 0, [U]. The VB PES has a spectral shape 
quite similar to that of the Ni 2^3/2 core- level spectrum. 
The only difference between the VB and the core-level 
spectra is the feature (which was enhanced in the 846 
eV spectrum). The character of the features A'-D' can 
be seen in more detail in Fig. 2(b). One can see that 
the features A' and B' are dominated by the 3d^Z_ state, 
with weaker contributions from the 3d^L state. The fea- 
ture G' is dominated by the 3d^ state. The feature has 
mainly 3d^L character. The most important point is that 
two lowest energy electron removal states mainly arise 
from the "bound state screening" 3d^Z_ state. We recall 
here that the standard single site GI model for the VB 
spectra assigns both the features A' and B' to the 3d^L 
state. Further, because the MSG is not yet reported for 
the VB PES of NiO, the present model alone reproduces 
both the core-level and VB spectra for a common set 
of parameters. Our results indicate that the feature A' 
and B' originate from the ZR doublet bound state. This 
indicates that the transport properties of the resistance 
switching observed in NiO and deficient NiO involves the 
ZR bound states, which has a composite character. This 
is because the first ionization state can be reached by 
simply removing one electron in VB PES. Two theoret- 
ical studies have identified the ZR feature for the first 
ionization state 13, lH. In fact, one study has general- 
ized it to the ZR bands for the first ionization state in 
GT insulating TM monoxides from MnO to GuO, as well 
as layered high-Tc cuprates[15]. 

The bound state screening mechanism is further con- 
firmed by resonant PES at the Ni L3 edge excitation. 
In Fig. 2 (c) we show the theoretical and experimental 
resonant PES for incident photon energies labeled a-d 
in the XAS shown in the inset. The off-resonant spec- 
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FIG. 3: (Color online) Schematic of the effect of final states 
for core-level 2p PES and VB PES 



found an improved agreement for the Ni 2p core- level, 
VB and resonant PES of NiO with a common parameter 
set. All features in the VB and the core- level PES were 
reproduced well by including the bound state screening 
effects, which is beyond the standard CI model. For NiO, 
the charge transfer effects from the ZR doublet bound 
states determine the first ionization state (i.e. the lowest 
hole-addition state). We believe that our results provide 
key features to understand the conduction mechanism 
of NiO and doped NiO. The present study provides the 
most complete picture to describe the core-level and the 
VB electronic structure of NiO, consistent with the spin- 
fermion model and the LDA+DMFT results. 

M. T. would like to thanks G. A. Sawatzky and A. Fu- 
jimori for useful discussions and illuminating comments. 



trum {hu=S4:6 eV) was discussed above in comparison 
with the hu=1.25 keV spectrum. Our on-resonant spec- 
tra are very consistent with a previous study [25]. The 
resonant spectra were calculated with (thick solid line) 
and without (dashed line) the bound state screening. The 
agreement between the theory including the bound state 
screening and the experiments is excellent for the energy 
dependence of the spectra. In particular, the spectral in- 
tensity changes of the two lowest energy features A' and 
B', corresponding to the ZR doublet x^—y^ and ?>z^—r'^ 
bound states, match nicely with the calculations. 

Finally, we discuss the similarity in the VB and the 
core- level PES. Figure 3 (middle panel) shows the Ni 
2p3/2 core-level HX PES and the 1.25 keV VB PES on 
the common energy scale. As is well known, the spectral 
shape presents common features [oj. In the final state of 
the core- level PES, a core-hole is created, and the energy 
of 2p^3(PL and 2p^3(PZ_ states are pulled down by the 
attractive core- hole potential —Udc with respect to that 
of 2p^3d^ states. Because A*<A, the 2p^3d^Z_ state lies 
below the 2p^3dPL state. In the VB PES, a 3d electron is 
removed in the final state. As shown in the upper panel 
of Fig. 3, the energy difference between final states are 
E{3d^L)-E{3d^)=A-U and E{3d^Z)-E{3d^)=A'^ -U , 
respectively. As a result, the 3d^L and the 3d^Z_ states 
are lower in energy than the 3d^ state. Again, the 3d^Z_ 
state lies below the 3d^L because A*<A. Therefore, the 
assignment of the feature A and A', as well as C and C^ 
are the same. But the assignment for the B and B' are 
not the same. The feature B originates in the 2p^3d^L 
while the feature B' originates in the 3d^Z_ and the fea- 
ture is the equivalent of the feature B (see Fig. 2(a)). 
This difference is understood in terms of larger multiplet 
splitting in VB PES leading to the bound state features 
A' and B', while the smaller multiplet splitting in core- 
level result in a narrow peak A. 

In conclusion, by using the extended CI model, we 
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